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CN ■ Abstract 

>: 

. We review the spontaneous magnetization of the vacuum of non-Abe- 

CN ■ Uan gauge fields at high temperature. The standard model of particles is 

investigated as a particular example. By using both analytic methods of 
quantum field theory and gauge field theory on the lattice, we determine 
the Abelian (chromo) magnetic fields in the restored phase of the model at 
CN ■ high temperatures T > Tew The fields are stable and temperature depen- 

dent, B = B{T). We investigate the mechanisms of the field stabilization 
in detail. The screening parameters for electric and magnetic fields - the 
^ ■ Debye, m£)(-B,T), and magnetic, mmagn.{^i'^)i masses - are calculated. It 

^ . is shown that, in the field presence, the former one is smaller than at zero 

field. The magnetic mass of the (chromo) magnetic fields is determined to 
be zero, as for usual U{1) magnetic field. We also show that the vacuum 
magnetization stops at temperatures below the electroweak phase tran- 
sition temperature, T < T^w ■, when a scalar condensate creates. These 
properties make reasonable a possibility that the intergalactic magnetic 
fields observed recently were spontaneously generated in the hot Universe 
due to vacuum polarization of non- Abelian gauge fields. We present a pro- 
cedure for estimating the field strengths B[T) at different temperatures. 
In particular, the value of B(Tew) ~ lO^^G, at T^^ is estimated with taking 
into consideration the observed intergalactic magnetic field Bq ~ 10~^^G. 
The magnetic field scale is also considered. Some model dependent pecu- 
liarities of the phenomena studied are briefiy discussed. 
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1 Introduction 



Strong magnetic fields present in all the objects of the observed Universe. As it 
is believed, these fields presented at previous stages of its evolution. They also 
are produced in the heavy ion collisions at high energies and influence various 
processes and characteristics of plasma, in particular, the deconfinement phase 
transition temperature. So, searching for effective mechanisms for generation of 
different type magnetic fields at high temperature is of great importance for either 
particle physics or cosmology. These mechanisms could serve as a theoretical basis 
for investigations of the QCD vacuum at high temperature in magnetic fields [1] , 
[2], [3], [5] and the primordial magnetic fields in the early Universe [6], [7], 
[8], [9]. Recent discovery of extragalactic magnetic fields possessing the field 
strength of the order Bq ~ IQ-^^G at 1 - lOMpc scale |10], [II] is one of the most 
interesting events of modern cosmology. In Ref. [I2] a model-independent 95% CL 
interval 1 • lO^^^G < -B < 3 ■ lO^^^G is determined, and the femtogauss values are 
actual field strengths in intergalactic space. First of all, it makes most reasonable 
the cosmological origin of primordial magnetic fields. From theoretical point of 
view, this discovery restricts in an essential way possible processes supplying the 
generation of fields in the hot Universe [H]. So, the searching for mechanisms of 
their creation is intensified. Of course, the prime candidates are the primordial 
fluctuations. But there are more candidates, for a review see [13], [H], [15]. A 
challenge is to produce coherent magnetic fields on different large scales in almost 
empty intergalactic space. In the present paper we review one of possibilities - 
spontaneous magnetization of the vacuum of non-Abelian gauge fields at finite 
temperature. This phenomenon is related with one of the basic properties of 
non-Abelian gauge fields - a large magnetic moment (gyromagnetic ratio 7 = 2) 
of charged gluons and M^-bosons. The fields generated spontaneously are long- 
range, temperature dependent and stable. All these features make them very 
attractive objects for various applications. 

The spontaneous magnetization of the vacuum of non-Abelian gauge fields at 
finite temperature was discovered recently in SU{2) gluodynamics. It was found 
by analytic methods in [16], [T7], [18], [19] and confirmed by lattice simulations 
in [20], [21]. The basic idea rests on the known fact that the spontaneous vacuum 
magnetization is the consequence of the spectrum of a color charged gluon, 

pl=pl + {2n + l)gB (n = -1, 0, 1, ...), (1) 

in a homogeneous magnetic background, B, described by the potential 

= Bx2 d.sS''^ (2) 

where a is weak isotopic index, and p\\ is a momentum component along the 
field direction. This Abehan type field is the solution to the classical gauge field 
equations without source terms. Here, a tachyon mode is present in the ground 
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state {n = —1). In fact, one observes < resulting from the interaction of 
the magnetic moment of the spin-1 charged particles with the magnetic field. 
This phenomenon was firstly discovered by Sawidy [22] at zero temperature, 
T = 0, and become known as Sawidy vacuum. However, at zero temperature, 
this state is not stable. It decays under emission of gluons until the magnetic 
field B disappears. 

This picture changes with increasing temperature where a stabilization sets in. 
The stabilization is due to a vacuum polarization. It depends on two dynamical 
parameters appearing at T 7^ 0. These are a magnetic mass of the color charged 
gluon, rrimagn., and a Ao-condensate, which is proportional to the Polyakov loop 
|23j . This field configuration is stable, its energy is below the perturbation one 
and the minimum is reached for the field of order gB ~ g^T^/ logT. This phe- 
nomenon is common for different SU{N) gauge fields which can be used to extend 
the standard {SU{2) x U{l))ew x SU{3)c model of elementary particles. 

An important property of such temperature dependent magnetic fields is the 
vanishing of their magnetic mass, rrimagn. =0. It was found both in one-loop 
analytic calculations |23] and lattice simulations [25]. The mass parameter de- 
scribes the inverse spatial scales of the transverse field components, similarly to 
the Debye mass related to the inverse space scale for the electric (Coulomb) 
component. The absence of the screening mass means that the spontaneously 
generated Abelian chromomagnetic fields are long-range at high temperature, as 
is common for the U{1) magnetic field. Hence, it is reasonable to believe that, 
in the hot Universe, at each stage of its evolution spontaneously created, strong, 
long-range magnetic fields of different types have been present. The fields influ- 
enced various processes and phase transitions. 

The aim of the present paper is to describe the main points of the vacuum 
magnetization phenomenon, investigate the properties of the created fields and 
discuss some possible physical consequences of it. In fact, this phenomenon is 
a non-perturbative one. So, the main results have been obtained both by the 
continuum field theory methods and by the lattice simulations, which occur in 
agreement with each other. Every method of the calculations mentioned has spe- 
cific features which will be described below. We consider the magnetization in 
the standard model (SM) as well as the minimal supersymmetric standard model 
(MSSM) at high temperatures and compare the results. We describe the mecha- 
nisms resulting in stabilization of the magnetized vacuum. As an application, we 
consider the origin of intergalactic magnetic fields which will be investigated in 
more details. We should note here that a difficult problem of magnetogenesis is 
to relate the field strengths generated in the early Universe with the present day 
fields, that depends on numerous factors and is model dependent. These points 
will be also discussed. 

In the next section, we start with analytic calculations and describe the con- 
sistent effective potential (EP) V(T, B, (p^) of magnetic and scalar fields used in 
studying of the spontaneous vacuum magnetization in the SM. The main obser- 
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vation here is that the magnetization does not happen when a sufficiently large 
scalar condensate 0c 7^ is present. It means, in particular, that after elec- 
troweak phase transition (EWPT) the spontaneous vacuum magnetization does 
not happen. That is why, we consider the magnetization in different models of 
particles for the case of 0c = 0, only. Sect. 3 is devoted to investigation of the 
vacuum magnetization at high temperatures for the SM. Sects. 4, 5 are devoted 
to calculation of the Debye and magnetic masses for neutral and charged gluons 
at high temperature in the external magnetic field presence. For this the one- 
loop gluon polarization tensors are calculated and investigated. Then we turn to 
non-perturbative methods of calculations and investigate the vacuum magneti- 
zation and the magnetic mass of neutral gluons at finite temperature in the field 
presence by applying Monte Carlo methods on the lattice in Sects. 6 and 7. We 
obtain the results coinciding with that of the analytic calculations. In Sect. 8 we 
develop a procedure to relate magnetic fields generated in the hot Universe with 
the present day magnetic fields and estimate the field strength at the EWPT 
temperature Tg^. In Sect. 9 the scale of magnetic fields at different temperatures 
is estimated. Discussion of the results obtained and prospects is given in the final 
section. 



2 Effective potential at finite temperature 

As we noted above, the spontaneous vacuum magnetization and the absence of 
the magnetic mass for the Abelian magnetic fields are non-perturbative effects 
to be determined, in particular, in lattice simulations [20], [25]. The main con- 
clusions of these investigations are that the stable magnetized vacuum does exist 
at high temperature and that the magnetic mass of the created field is zero. 
Concerning the actual value of the field strength, it is close to the one calcu- 
lated within the consistent effective potential which takes into account one-loop 
plus daisy diagrams. So, in analytic calculations we restrict ourselves with this 
approximation. 

The complete EP for the SM is given in the review [26j. Its general structure 
reads 

V{T, 5, (t>c) = ^ + V'''\<\>c) + V^^\T, B, 0c) + y""5(r, B, 0c), (3) 

where 0c is a scalar field condensate. The ffist two terms describe the tree level 
contributions of classical magnetic and scalar fields, V^^\T, B, (pc) is the one- 
loop contributions of all the fields and the last term presents the contributions of 
ring (or daisy) diagrams giving the main long-range correlation corrections. To 
find V^^^^{T, B,(j)c) the one-loop polarization operators of charged, Il{T,B) and 
neutral, 11° (T, B) gluons in the external field and at finite temperature have to be 
calculated in the limit of zero momenta. Then V"^^{T, B) is given by the series 
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Figure 1: The neutral gauge fields and scalar field ring diagrams giving contri- 
butions to the effective potential. 




Figure 2: The charged gauge field ring diagrams giving contributions to the 
effective potential. 

depicted in figures [Tl[2l Here, the dashed lines describe the neutral gluons and the 
wavy lines represent the charged ones, blobs stand for the one-loop polarization 
operators. The diagrams with one blob show the two-loop terms of the EP, with 
two blobs - three-loop ones, etc. 

The one-loop terms have the orders aem = e^/(47r) and A in a fine structure 
constant and a scalar self-interaction coupling. The order of the rings is Og^, A^/^, 
correspondingly. So, these diagrams are dominant compared to the two-loop 
contributions having the order a^^, A^. As it is well known, the EP ([3]) is a 
consistent approximation because it is real at sufficiently high temperature (see 
for details [26j and references therein). So, it will be used in the following. 

In the present study we are interested in two limits of it, 

1. weak magnetic field and large scalar field condensate, eB < (pl.T ^ 0, 

2. the case of the restored symmetry, 0c = 0, i? 7^ 0, T 7^ 0. 

For the former case we show the absence of spontaneous vacuum magnetization 
at finite temperature. For the latter one we estimate the field strength at high 
temperature. 

To demonstrate the first property we consider the one-loop contribution of 
ly-bosons calculated in a standard way as the sum over the spectrum ([1]) (see 
Ref. [26] for details): 

Vi'\T.K^) = :^,t f'~f''^ K,(nP{<p--hfl^) (4) 

Th 
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n labels discrete energy values and Ki{z) is the MacDonald function. Here and 
in what follows we use the dimensionless variables: h = eB/M^, = 0/0o (0o is 
the value of scalar field condensate at zero temperature), and /3 = M^/T, M^, is 
the H^-boson mass. In fact, other terms of the EP have to be added. However, 
to elucidate the role of the scalar field condensate this term is sufficient. 

Now, let us show that the spontaneous vacuum magnetization does not happen 
at finite temperature and for non-small values of the scalar field condensate 7^ 0. 
To do that we notice that the magnetization is produced by the gauge field 
contribution, given in Eq. (jl]). So, we consider the limit of |§ <C 1 and (jP' > h. 
For this case we use the asymptotic expansion of Ki{z), 

K,{z) ~ .f^e-^ (1 + -- + ■■■), (5) 

where z = n(3{(j? ± hY^. Now, we investigate the limit of /3 — ?■ 00, <C 1. For 
this case the leading contribution is given by the first term of the temperature sum 
in Eq. (jl]). We can also substitute {cj^rLhY/'^ = 0(1 ±2^). In this approximation, 
the sum of the tree level energy and (j4]) reads 

2 n^/" ^'"4,"' \ 234,) ' 

The second term is exponentially small and the stationary equation = 
has the trivial solution /i = 0. This estimate can be verified easily in numeric 
calculation for the total EP. Hence, we conclude that after symmetry breaking the 
spontaneous vacuum magnetization does not take place, as at zero temperature 



The main goal of our investigation is the restored phase of the SM. So, we 
adduce the high temperature contribution of the complete effective potential 
relevant for this case using the results in Ref. |26j. Firstly we write down the 
one-loop ly-boson contribution as the sum of the pure Yang-Mills weak-isospin 
part {B = 5(3)), 



where g is weak isospin gauge couphng constant, /x is a temperature normalization 
point, and the charged scalars [T8] . 



V}PiB,T) = + ^^-^^^ + Oig^B% (8) 

4o TT^ /i iZ TT 

describing the contribution of longitudinal vector components. Remind that elec- 
tric charge is expressed through g and the Weinberg angle as e = gsinOw The 
first term in Eq. ([7]) is the tree-level energy of the field. This representation 
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is convenient for the case of extended models including other gauge and scalar 
fields. Depending on the specific case, one can take into consideration the parts 
(Ej), dHl), correspondingly. In the SM, the contribution of Eq. ([H]) has to be taken 
with the factor 2, due to two charged scalar fields entering the scalar doublet of 
the model. In the case of the Two-Higgs-Doublet SM, this factor must be 4, etc. 
The imaginary part is generated because of the unstable mode in the spectrum 
([1]). It is canceled by the term appearing in the contribution of the ring (daisy) 
diagrams for the unstable mode p!9]. 



gBT r ~ ~ii/2 XgBf/^T 



unstable 



2tt 



UiB,T,n = -l)-gB + . (9) 



Here 11(5, T,n = —1) is the mean value for the charged gluon polarization tensor 
taken in the ground state n = — 1 of the spectrum ([1]). If this value is sufficiently 
large, spectrum stabilization due to radiation correction is taking place. This 
possibility formally follows from the temperature and field dependence of the 
polarization tensor in the high temperature limit T — )■ oo [2S]: Il{B,T,n = 

— 1) = c g'^T\J gB, where c > is a constant which must be calculated explicitly. 
At high temperature the first term can be larger than gB. Hence we see that 
the accounting for of rings leads to the vacuum stabilization. It is obvious in the 
high temperature limit. But it also is the case for non- asymptotic temperatures. 
Really, in the latter case, the imaginary part will be smaller than in Eq. ([7]). 
The high temperature limit of the fermion contribution looks as follows [26] . 

a 1 ~ T 

Vf ermion = J2 (^0) 

TT ^ 6 /i 

where the sum is extended to all leptons and quarks, and qj is the fermion electric 
charge in positron units, /i is a normalization point. Hence it follows that in the 
restored phase for fi ^ T^w all the fermions give the same contribution. 

Now, let us present the EP for ring diagrams describing the long-range corre- 
lation corrections at finite temperature [29], [30] . 



, (noo(o)) 

327r2 



( 471 \ 3_ 
H/3(noo(0))V2;+4 ^. 



where the trace means summation over all the contributing states, Hqo = H0(fc = 
0,T,B) for the Higgs particle; m|, = Hqo = ^ooik = 0,T,B) are the zero-zero 
components of the polarization functions of gauge fields in the magnetic field 
taken at zero momenta, called the Debye mass squared, 7 is Euler gamma. These 
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terms are of the order ~ 5'^(A^/^) in the couphng constants. The detailed calcula- 
tion of these functions is given in Ref. [26]. We give the results for completeness, 



n<^(o) 




6e2 



+ 



3e3 



sin^(2^VK 

^ - \qfB\K^ 



sin^ 9w 



(12) 



+ 



{eB) 



1/2 



Svrsin OwP 



3y2C(- 



1 h 

'2' 2' 



Here Kj = = '^^^2°""' and A is the scalar field coupling. The terms 
give standard contributions to temperature mass squared coming from the bo- 
son and fermion sectors. The i?-dependent terms are negative (the value of 
(3a/2(^(— |, |) = —0.39). They decrease the value of the screening mass at high 
temperature. The Debye masses squared for the photons, Z-bosons and neutral 
current contributions are, correspondingly. 



m 



D,7 



m 



D,Z 
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^ ' tan^ Ow + 



+ 0{eB(3^ 
1 



3/32 ' " 



(13) 



4cos2^vk/ 3/32 



+ 0{eBp'), 



m 



9' 



D , neutral 



8 cos2 OwP"^ 



(l + 4sin^^VF) +0(eS/32). 



As one can see, the dependence on B appears at order 0{T 2). 

The py-boson contributions to the Debye mass of the photons are 



^'dw = ^9^ sin^ Ow 



3/32 



(g sinOwB) 
27r/3 



1/2^ 



(14) 



Interesting feature of this expression is the negative sign of the next-to-leading 
terms dependent on the field strength. Finally, we give the contribution of the 
high temperature part in Eq. (Q Il{B,T,n = —1), 



U{B,T,n = -1] 



-a 



3.26 {g sine wBy/^ .{g_smewBfl'^ 



(15) 



This expression was calculated from the one-loop VT-boson polarization tensor 
in the external field at high temperature (see Eqs. (|7Tll . (1701) below). It is one 
of the main objects for the problem of interest. So, we consider the calculation 
of both neutral and charged modes of it in Sects. 4, 5. The expression f|T5|) con- 
tains the imaginary part which comes from the unstable mode in the spectrum 
(II]). Its value is of the order of the usual dumping constants in plasma at high 
temperature. It will be thus ignored in actual calculations in what follows. In 
fact, this part must be calculated in a more consistent scheme which starts with 
a regularized stable spectrum. The stabilization can be realized not only by the 
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radiation corrections but also some other mechanisms. We observed a stable vac- 
uum state in the lattice simulations [20]. Therefore, we believe that this problem 
has a positive solution. Other important observation obtained is the absence 
of the vacuum magnetization at not small values of the scalar field condensate, 
(pc 7^ 0. In particular, this means that just after the electroweak phase transition 
this phenomenon does not hold. 



3 Vacuum magnetization at high temperature 

The EP (j4]) (and other terms in Ref. [26]) is expressed through the well known 
special functions. Therefore, it can easily be investigated numerically for any 
range of parameters entering. As usually, it is convenient to introduce the di- 
mensionless variables: the field $ = {gBf/'^/T and the EP g) = V{B, T) /T^. 
The vacuum magnetization at high temperatures, T >> (gBY^'^, $ — )■ 0, can be 
investigated within the following limiting form of the total EP, 

where ~ g'^ terms were omitted. The logarithmic term is signaling asymptotic 
freedom in the field at high temperatures. In this approximation, the role of 
the ring diagrams is reduced to the cancelation of the imaginary part in Eq. (jH]), 
only. The first term has the order ~ g^^* and can be neglected. This follows from 

the high temperature estimate for n(T, n = — 1)|t->oo ~ 9^T\fgE (see Eq. ( 17T1) 
below). Thus, in given approximation the EP is real. From Eq. ( 1T6|) we obtain 
for the minimum position 

2^ i + (|-E/|g^)^iog| 

Hence, we come to the conclusion that the ferromagnetic vacuum state does exist 
at high temperatures. The field strength flTTl) is proportional to the coupling 
constant g^. Such type dependence is opposite to the case of the usual sponta- 
neous breaking of symmetry where the value of the scalar condensate depends 
non-analytically on the coupling: 0o ~ V VX. 

Let us discuss in details the stability of the condensed field. First let us note, 
for the specific value of the field strength ( ITTl) the order of the terms in Eq. (fTSll 
is ~ g^, that approves the consistency of the approximation used. Further, the 
second derivative of the EP is positive for gBc that means we have a minimum. 
The field is not changing in the direction a = 3 of the isotopic space. To check 
is this the case or not for the perpendicular directions = 1, 2 or a = which is 
responsible for excitation of charged fields W^, one has to calculate the effective 
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mass squared M'^{Bc,T). Here, we consider the one-loop case. Substituting the 
value ^S(i) Eq. ([HD in the one- loop polarization function Eq. ( !T5|) . we find that 
the effective mass squared, M'^{B^^\T) = -[ReU^^\Bi^\ n = -1) - gB^^^] > 0, 
is positive for any negative Rell at sufficiently high temperature. This is because 
the second term is logarithmic suppressed (~ l/log(T/yu)) compared to the first 
one. Thus, the vacuum stabilization is expected in this consistent calculation. 

The above analytic investigation unambiguously determined the possibility of 
the vacuum magnetization at high temperature, although a number of questions 
has to be studied in order to derive a final picture. In fact, the radiation correc- 
tions are not the only candidate for stabilization of the vacuum. Other factor is 
a so-called Ag-condensate which is proportional to Polyakov's loop. This param- 
eter also acts as a stabilizing mass term at high temperature [16], [31], [21]. Its 
role will be discussed in what follows. Let us note once again that we are dealing 
with infrared properties of gauge fields at high temperature and non-perturbative 
methods of calculations have to be used. So, below we investigate the magne- 
tization as well as the magnetic field parameters by means of analytic methods 
and simulations on the lattice. 

4 Temperature masses in the field presence 

As it is well known, at finite temperature particles acquire the temperature de- 
pendent masses which modifies interactions in plasma [321 [33]. At zero fields, 
this mass is of the order ~ gT (m|, = n44(T, ^4 = 0, /c — >■ 0)) for quarks and 
longitudinal gauge field excitations, where g is gauge coupling. For transversal 
gluons, this mass is of the order ~ g^T = Ilij{T,k4 = 0,k 0)). So, 

both chromoelectric and chromomagnetic fields become screened in hot QCD. 
These basic characteristic parameters of quark-gluon plasma (QGP) are very 
important. They determine the screening properties and plasma excitations of 
QGP for different type gauge fields and have been calculated in analytic quan- 
tum field theory at finite temperature and simulations on the lattice. Note that 
in scalar and spinor QED the longitudinal (electric, Debye) mass has the order 
~ eT and the magnetic mass of photons equals to zero [32]. Therefore the latter 
field remains a not screened long range one as in the vacuum. 

At non-zero external magnetic fields, these characteristics are not well investi- 
gated at present. As it was discovered in analytic calculations [M] and confirmed 
in simulations on the lattice [25] the magnetic mass of neutral gluons calculated 
in the external chromomagnetic Abelian field B = const described by the poten- 
tial A^^ = B xi 5"'^5ij,2 equals to zero. So that this field behaves as usual magnetic 
field in QED. It is long range one, in contrast to the zero field case. The Debye 
mass of these fields is field-dependent, and its value decreases with growth of the 
field strength [23]. So, the color Coulomb field becomes a more long range one 
compared to the zero field case. 
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The Debye mass of charged gluons (calculated in SU{2) gluodynamics) is also 
field-depended and coinciding with the neutral gluon value [34J- As concerns the 
magnetic mass for charged gluons, it remains not well calculated. Only rough 
estimates have been obtained already [19]. But this parameter is of paramount 
importance for plasma. 

From dimensional analysis it follows that the PT at high temperature has 
the form ng^ (5,T,p3 0) ~ c^/gBT, where c is a constant which has to be 
computed. So that if c < the spectrum stabilization at high temperature 
(T — )■ oo) holds. If c > 0, the spectrum and, therefore, external field remains 
unstable as at zero temperature. 

Hence one also can check a stabilization mechanism of the spectrum due to 
radiation corrections in the field presence. If the last term occurs to be large 
positive value, the spectrum stabilization takes place at a given temperature. 
Even a general sign of the last term is essential. This is because in non-Abelian 
gauge theories the sign of the radiation corrections usually is not changed when 

2 

next in the a = parameter corrections are taken into consideration. Therefore 
a possibility for stabilization will be discovered even at a one-loop level. 

In fact, we have to distinguish two situations: 1) the spectrum at high temper- 
ature in the external field presence which can be positive for any field strengths 
at high temperature for c < 0; 2) the spectrum in the field which is spontaneously 
generated at high temperature {B = B(T)). In the later case the value of the 
gluon effective mass 

mlff^ = -gB -{n = -l|n^f (^4 = 0,5,r,p3 ^ 0)|n = -1) (18) 

can be negative dependently on the value of B(T). The fist case is of interest, 
for example, for the QGP in pp-collisions, when one can expect strong external 
chromomagnetic fields at high temperature produced by quark currents. The 
second case is of interest for cosmology at T > T^. 

5 Polarization tensor of neutral gauge fields at 
high temperature 

In this section, we investigate the properties of the neutral gauge fields in the 
restored phase of the SM by analytic methods. As before, in actual calculations 
we concentrate on the electroweak sector SU{2) x f/(l)y of gauge fields. The 
color SU{?))c gluons posses similar properties which, for brevity, we note, only. 

First problem to study is space scale of gauge fields at finite temperature 
in the background field presence. To investigate this problem, the polarization 
tensor of neutral gauge fields n(T, 5, ^4, /c), where ^4, /c are gluon momentum 
components, must be calculated. Before doing that, we remind the structure of 
the gauge fields in the SM. This is necessary, because at high temperature the 
SU{2){SU{?>)c) neutral components of gauge fields are generated spontaneously. 
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Figure 3: Basic graph for polarization tensor (left figure) and graph with one 
vertex and a closed line (right figure). 



The polarization tensor consists basically of the graphs shown in Fig. |3] where 
gluon and ghost lines must be inserted. 

The electromagnetic and Z-boson potentials read 

= TPTF^^'^^^^^'^^' ^^^^ 

^ 3 



where A^, are the Yang-Mills gauge field third projection in the weak isospin 
space and the potential of the hypercharge gauge field, and g and g' are SU{2) 
and U{1)y couplings, correspondingly. At high temperature, 7^ 0, the field is 
spontaneously generated. The hypermagnetic field is not spontaneously gener- 
ated, = 0. After the electroweak phase transition, the Z-boson acquires mass 
and the field is screened. The only component = —^==g'Aj^ = sin 9^^^^ 

remains non-zero. Here 9]y is the Weinberg angle, tanOw = ^. Hence, at low 

temperature the remnant of the field B is the magnetic field B = sin OwB. Be- 
low, for brevity, we write B instead B for the background field. The difference 
will be taken into account when necessary. 

As it is well known, at finite temperature particles acquire the temperature 
dependent masses which modifies interactions in plasma [521 • At zero fields, 
this mass is determined as mjy = H44(T, ^4 = 0, — )■ 0) for quarks and longitudi- 
nal gauge field excitations, and has the order rriD ^ gT. For transversal gluons, 
this mass is of the order rrimagn. ~ g^T {m^^g^ = {s\Ilij(T, = 0, A; — )■ 0)|s), 
where \s) is the transversal state of the gluon). So, both chromoelectric and 
chromomagnetic fields become screened in hot SM. These basic characteristic pa- 
rameters of gauge fields have been calculated in analytic quantum field theory 
at finite temperature and simulations on the lattice. Note that in scalar and 
spinor QED the longitudinal (electric, Debye) mass has the order ~ eT and the 
magnetic mass of photons equals to zero [32]. Therefore the latter field remains 
a not screened long-range one as in the vacuum. 

This picture is modified in the background field presence. In this case it 
is reasonable to consider color gluon field, (a = 1, 2, 3), separately in terms of 
charged W^-bosons (gluons), Wj^ = -^{V^ ± ^^1^), and neutral, = A^^, ones. 
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Figure 4: The neutral polarization tensor. 

The one-loop neutral gluon polarization tensor (PT) in the background field 
at finite temperature was calculated and investigated in Ref. [21]. An imaginary 
time formalism was used and, in particular, both the parameters, mij{B,T) and 
'^magn.iTi B) havc becu derived. 

In general, the PT can be present in the form [24] : 



10 

n,v(A:) = E ^^Hk,B,T)T^l, , (20) 

i=l 

where T^}^, are ten structures out of the momentum fc^, medium velocity and 
5x\i and the form factors 11'^*) (fc) depend on the external momentum through 
the variables = k\ + fc| and h'^ = kf + A;| at zero temperature and /i^, k^ and 
fcs at finite temperature. These structures form a complete tensor basis. 

The one-loop PT has the following representation in momentum space (see 

Fig.SD 

Uxx'{k)=T Yl J T^Il{p,P4,k,k^)xx', (21) 



Af=-oo_ 



where in the integrand we noted explicitly the dependence on the external mo- 
mentum and the momentum inside loops. The integrand looks as follows 

Ilxy{p,k) = T^yxG^^\p)T^i^:x'Gu'u{p- k) (22) 
-pxG{p){p - k)x'G{p - /c) - (p - k)xG{p)px'G{p - k) 
+Gxx'{p) + Gx'x{p) - 25AA'TrG'(p). 

The second line gives the contribution from the ghost loops and the third one is 
due to the tadpole diagrams. The vertex factor is 

r^^^A = ^^.v{k - 2p)x + 6x^^{p + k)^ + 6xu{p - 2k)^. (23) 

It should be remarked that all graphs and combinatorial factors are exactly the 
same as in the well known case without magnetic field. On this level the only 
difference is in the meaning of the momentum = idf^ + iB^^ which in our 
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case depends on the background magnetic field. The four particle vertexes are 
momentum independent and have the same form as at zero external field. 
The propagators are given by 



1 2 
G{p) = -= ds e-'^ , 
V Jo 



{p - kf 



G{p -k)= . = r' dt e-*(P-'=)' (24) 



for the scalar lines and by 



1 

p2 ~2iF 



Gxx'ip) = , , =/ dse-'^E-^ 



XX' 



OO 2 



AA'; 



|2 + 2iFy^^, io 



Gxx'{p-k) = . TTrr^] = dt e~'^^-^^ E-^i, (25) 



{p-kf + 2iF^ 
for the vector lines with 

^AA' = (e'^'^),,, = 4' + ^^AA' sinh(2.) + cosh(2s). (26) 

Here and below, for simplicity, we set the field strength gB = 1. That means 
we measure all quantities in units of gB. To return to the dimensional variables 
one has to substitute s gBs, etc. 

At zero temperature, the momentum integration can be carried out by means 
of Schwinger's algebraic procedure [35] and converted into an integration over 
two scalar proper-time parameters, s and t. Here we adduce the known results 
in order to present their modifications at T 7^ [21]. The basic exponential is 

e = e-'P'e-'^P-''^' (27) 

and the integration over the momentum p is denoted by the average (...). The 
following formula holds: 



exp 



-k{jh^^^ + ^S-)k 



(47r)2(s + t)sinh(s + t) '''^^^ 

with S = tanh(s) and T = tanh(i:). 

Within the developed formalism, the polarization tensor becomes the expres- 
sion of the type 

/•OO /"OO 

Uxyik)= / ds dt Mxx'ip,k){eTis,t)), (29) 
JO Jo 

where in Mxx' [p, k) we collected all factors appearing from the vertexes and from 
the lines except for that going into 0T(s,t) (Eq. (59) in Ref. [2] 



+°° / JV^ k tN \ 

{e{s,t))T = E exp - ^^^^ , ,, +'T^hvF] (30) 



Af=-oo 
+00 

E 0T(^,t). 



4T2(s + t) (s + t)T^ 



N=-oo 
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We introduced the notation 6^(5,^) which is the basic function appearing in 
all form factors. The function (G(s,t)) is given in Eq. (!28|) . In this way the 
finite and zero temperature formalisms are related. All the necessary formulas 
for momentum integration are given in Ref. p^ . 

Now, we calculate the magnetic mass of neutral gluons. It can be determined 
in the imaginary time formalism through the PT, 

^liagn.is) = (s | n(fc, 5, T) | s)fc,=o,fe2^0 • (31) 

Here, = 2tcNT, N = 0, ±1,±2, ... is a Matsubara frequency, k"^ = P + 
is a momentum squared. The mean value is calculated in the two states of 
polarization (denoted as s = 1 and s = 2 [36], [21]) transverse to the gluon 
momentum A;^. 

The factors M*^*-* giving contributions to the matrix elements of 11 for the 
infrared limit of interest in Eq. (!3T]) to the states s = 1, s = 2 are [23]: 

. 1 — cosh(o) cosh(f ) „ ^ , , , , 
^2 = 4 ^l^^-^-2 + 8cosh(g)cosh(0, 

Ma = -2cosh(2g) ^^^^|^j^| - 2 + 6cosh(n cosh(g), 
gsmh(g) 

Ms = -2 + 2cosh(g)cosh(0, (32) 

where q = s + t, ^ = s — t = q{2u — l) and s = qu, t = q{l—u). The contributions 
to the particular polarization states s are (Eq. (142) Ref. [24j): 

{s = l\U{k)\s = 1) = h^U2, 

(s = 2|n(A;)|s = 2) = /^^(ng + ns). (33) 

We have to calculate the form factors 112, II3 and lis. 

These matrix elements are the product of h'^ and expressions which have a 
finite limit for h"^ = 0, 

n, = nf + 0(/i2). (34) 

The quantities nf'^ can be calculated analytically in terms of simple zeta func- 
tions. 

From the above expressions (!29|) - (!32|) . in leading order for T — )■ 00, which 
picks just the = 0-contribution, we note 

In these expressions, the integration over u can be carried out exphcitly, 

- L- M-(a) (36) 

~ (47r)3/2y^y sinh(g) ^^^^^^^ ^^^^ 
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with 



M2(g) 
MM 
M,{q) 



4 



coth(g) + . , . .2 



4 4 , , , 
H — sinn(2g), 



2 3 

-2 cosh(2g) (-1 + g coth(g)) + - sinh(2g), 

q 

-2 + -sinh(2g), 
9 



(37) 



and the g-integrations remain. 

These expressions are formally divergent for g — )■ oo. This divergence results 
from the tachyonic mode. Here we have to remember that all formulas above are 
written in Euclidean representation (basically, for technical reasons). In fact, we 
have to start from the Minkowski space representation which can be reached by 
an 'Anti'-Wick rotation, g — > qe™!"^ . In the Minkowski space representation the 
parametric integrals are convergent using the usual 'ie'- prescription. But then, 
the contribution from the tachyonic mode in the loop cannot be Wick-rotated 
since, in momentum space, the corresponding pole is on the 'wrong' side of the 
imaginary axis of the momentum po- However, it can be 'Anti'-Wick rotated 
delivering a exponentially fast converging integral. The remaining part can be 
Wick rotated as usual. In this way, if starting from the Euclidean representation, 
the tachyonic part must be 'Anti'-Wick rotated twice, g — ?• ge*'^. The remaining 
part can be kept as it is. The subdivision into tachyonic and remaining parts 
must be done according to the behavior for g — t- oo. There is a freedom left of 
redistribution power like contributions. It can be used to avoid singularities in 
g = 0. 

We make the following subdivision into tachyonic parts (parts A), 



Mi 



Mi 



sinh(g). 



^ 2(1 -g) ^ 3^ 







sinh(g). 



sinh(g). 



(38) 



and remaining parts (parts B), = Mi — Mf". 

Being inserted into (!36l) . after the 'Anti'-Wick rotation, the A-parts constitute 
simple integrals which can be done immediately (for the moment we drop the 
common prefactors). 



n(o),A 


= 4iA/7r, 


Hf-^ 


= SiA/vr, 


n(o),A 


= iy/n . 



(39) 
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The integrals in the B-parts are directly well convergent. Their calculation is a 
bit more difficult, but after a number of transformations all can be taken into a 
form to be found in tables. The results are 



(0),B _ , , V" ' -v-'V-'V2 



3-3y2-87r + 2y27r)C(r 
in + 



6-6v/2 - 47r + 2 V2vr C(| 
njj = 5v7r + ^ 



n{o)-B 



^v^(2+(-4 + v^)C(^)) . (40) 
Together with the A-parts we get the corresponding numerical values [25]. pi], 

;-5.80 + 7.09i), 

;i.04-8.9i), 

;-4.21 + l.Sz) . (41) 



nf = 




T 




s/gB 


= 




T 


(47r)3/2 , 


/gB 


= 


g' 


T 


(47r)3/2 , 


/gB 



The above expressions have to be used in Eq. fl33|) to obtain final result. 
The sum of II3 + II5 equals, lis + lis = [—3.17 — 7.09z]. The imaginary part 
is signaling the instability of the state because of the tachyon mode, and the 
real one is responsible for the screening of transverse gluon fields. The real and 
imaginary parts are of the same order of magnitude. This is similar to the case 
of Landau's damping at finite temperature. 

Let us turn to the real part and substitute it in the operator Schwinger-Dyson 
equation 

D-i(P) = A;2 -n(fc) (42) 
for the neutral gluon Green function. We obtain for the mean values 

{s = l \D-\h^)\ s = l) = h^- Re(n2) (43) 

T \ 



1 + 5. 



/gB^ 



and 



s = 2 \D-\h^)\ s = 2 ) = -ReiUs + U^) (44) 



h'^(l + 7.09 ^ 



gB. 



These are the expressions of interest. 
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Two important conclusions follow from Eqs. f H3|) .f H4|) . First, for the trans- 
verse modes in the field presence, there is no fictitious pole similar to that of in the 
one-loop approximation for zero external field background at finite temperature 
[32] . [33], [37]. The external field acts as some kind resummation removing this 
singularity. Second, there is no magnetic screening mass in one-loop order. The 
transverse components of the gluon field remain long-range in this approximation, 
as at zero external field [32] . 

Possible resolutions of the zero one- loop magnetic mass are obvious: 1) the 
mass is generated in some kind resummation of perturbation series (as this is 
well known at zero external field case [38]); 2) there are no magnetic mass for 
neutral gluons as in the case of usual magnetic fields. The problem requires 
non-perturbative methods of computation. 

Second characteristic to calculate is the Debye mass determined in a standard 



way: m|) 



D 



-n44(T,5, h 
1 



47r2 



= 0,k^0) [Sg. It reads 
°° Bcos\i{2Bq) 



dq ^ 



sinh(i?g) 



3 



B 

4T2 



(45) 



where the dimensional parameters are restored. The function /(s) is dimension- 
less and it depends on the dimensionless variable s = It is chosen to satisfy 
/(O) = 1 such that it describes just the change which comes in from the magnetic 
field. 

The function f{s) can be easily computed numerically. Also, its asymptotic 
expansion is easy to obtain. We consider small s, i.e., high temperature, and 
represent /(s) in the form 



f{s) 



6 



1 cosh(2g) 
q sinh(g) 



(46) 



Here the first term in the square brackets delivers the zero field limit, /(O) = 1. 

The detailed calculation of it are carried out in Ref. [21], for the function f{s) 
we obtain 



f{s) 



1 + 



TT 



V2 - 1 



Hence, for the Debye mass we derive 

2 



D 



3^ 



-i 0.4775 



0.8859 



+ 



2T 



s + 



25 03) 

4 71^ ' 



2tt 




+ 0.4775 ( 



(47) 



(4J 



where the general factor g"^ and numeric values of the coefficients in Eq. ( H7|) were 
substituted. 
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In this way, from Eq. (1461) the Debye mass can be determined in terms of 
the Riemann Zeta-function. As expected, it has an imaginary part which results 
from the tachyonic instabihty. We note that the Debye mass in the magnetic 
field is smaller than without it. That is, in the field presence a color electric 
field becomes a more long-range one compared to the zero field case. From 
Eq. f HHj) it follows also that the imaginary part is next-to-leading order compared 
to the real one. The ratio p = Im m|)(i?)/Re m'jj{B) ^ 1. This property, in 
particular, means that that perturbative methods of calculations are sufficient for 
determining this parameter. In contrast, for the magnetic mass case the stability 
parameter p = Im(n(i?))/Re(n(i?)) is of the order ~ 1. In the latter case, p has 
the order of usual Landau's damping parameter in plasma. 

Similar behavior is proper to the neutral gluon field in the SU{3)c case. 



6 Polarization tensor of charged gauge fields at 
finite temperature 

Let us investigate the PT of charged gluons (VT-bosons) in SU (2) gluodynamics 
in the background field ([2]) at finite temperature [28], [39]. In the field, the 
polarization tensor can be constructed out of the vectors and d^, 



(49) 



where the third vector is rf^ = F^j^^p^, and we note px = l\ + hx, and the matrixes 



( ° ^ 




( Pi \ 




( 


P2 


\ 





, K = 




, d^^ = 




-Pl 




















V P4 J 




I ) 












/ 



10 

V 1 y 



■^fiX 



/ 1 \ 

10 



V / 



fiX 



f 1 \ 
-10 


V 0/ 



■ (50) 



Here, px is a momentum in the external magnetic field. It is an operator and it 
obeys the known commutation relation 



[Pi^,Pv 



(51) 



There are no more vectors or constant matrices, the polarization tensor may 
depend on. For instance, this reflects the residual rotational symmetry in the 
(p3,P4)-subspace. Hence the operator structures can also be constructed out of 
these quantities only. We mention that with our choice -8 = 1, F^^, in fl50|) is the 
field strength of the background field. 

Before writing down the decomposition of the PT in terms of form factors we 
mention one property which follows directly from the basic commutator relation 
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p-k 



Figure 5: The charged polarization tensor. 



the momentum obeys, namely = {p'^Sw + 2iFxy)py. As a consequence, 
for a function of the relations 

f{p>x = px'fip' + 2iF)yx (52) 

hold where now / must be viewed as a function of a matrix so that it itself 
becomes a matrix carrying the indices A and A'. The same is true with h"^ in 
place of p"^. 

The decomposition of the polarization tensor can be written in the form (see 

Fig. ED 

Uxx'ip) = E ^«(/^ + 2iF)xx" T§y + n°T°,. (53) 

i 

In general, the sum includes ten structures T^y, and T^x' = uxu'x defined in 
Refs. [36], [28]. The form factors U'''\f,h'^ + 2iF)xx> depend on f and only 
(besides their dependence on the matrices in fl50|) ). In fl53l) the form factors can 
be placed also on the right from the operator structures applying both relations 



The calculations of the charged PT differ in an essential way from the neutral 
case although all the steps and methods are similar. The main difference consists 
in the dependence of form factors on the operators h"^ + 2iF . We use the repre- 
sentation of the polarization tensor given in [28]. In momentum representation, 
the initial expression reads 

/dk ( 
J^^:^ YXupGu,y'{p-k)Tx'u'p'Gpp'{k) 

+ip - k)xG{p - k)kx'G{k) + kxG{p -k){p- k)x'G{k) 

+nl\lP°' , (54) 
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where the second hne results from the ghost contribution and the tadpole contri- 
bution is given by 



n 



tadpol 
AA' 



dk 

(2^ 
dp 



5xx'G,,{k) - Gik)xx 
Sxx'Gppip) + Gx'xip) - 2Gxx'{p) }■ 



(55) 



The vertex factor and propagators are given in Eqs. ( l23i) - ( 125|) . 

To proceed further, it is convenient to turn the momentum and the gluon field 
in color space into the charged basis, 



= Pa with 5 



I 



\ 



1 



1 1 

i —i 






1 
1 



(56) 



fj,a 



we have a color neutral and a color charged field, both of spin 1. The charged field 
will occupy in the presence of the background field Landau's levels and we expand 
with respect to the corresponding eigenf unctions. In this space, the momentum 
Pp of the charged gluon is an operator, whose components fulfill the commutation 
relation 

/ -1 \ 



[Pa,Pl3] 



iFn,R = i 



10 
10 
V 1 / 



B. 



(57) 



0/3 



As before, in this expression we write B instead gB. In what follows, we also put 
S = 1, for short. 

A basis in this space is given by the vectors | n,a)p_, Eq. (30) in [36j. The tree 
level energies of these states are 



En = ll + ll + B{2n + 1 + 2a) 



(n 



0,1 



a 



±11 



(5^ 



where I3 and U are the momenta in parallel to the background field and imaginary 
time, respectively. The lower state is n = with a = —1. It is tachyonic one. In 
this section, in contrast to Eq. ([T]), we use other system for marking the charged 
gluon states in the field. In the charged basis, the lower {tachyonic) state is 



t)a =\ 0,-1), 






V / 



0), 



(59) 
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where | 0) is the lower Landau level which is annihilated by the operator a in 



Pa 



—la 
h 

V h ) 



(60) 



and we note 

vA^)^ = vi\t)<. = ^. (61) 

Below, we will also use the notations = /| + l\ and h? = j)\ = aa^ + a^a. 

Now, we continue with using the representation of the PT as given by Eq. (51) 
in [28]. It results from the proper time representation of the propagators and the 
vertex factor given in Eqs. fl2^ . ( |25|) . ( l23l) and integration over k in Eq. ( 15^ . 

The one-loop PT in the imaginary time formalism can be written in the form 
of the integral over two proper time parameters and temperature sum 



Hav = E / ^^^^ 



tadpol 



(62) 



with M^y coming from the main gluon one-loop diagram, Mf^/ is the correspond- 
ing contribution from the ghost loop. The last term presents the contributions 
of the tadpole diagrams. Here we introduced the basic average 



4(s + t)T2 
exp(— if) 



+ 2s(up) e(s,t) 



eT(/^ h'^) = exp I 

where 0(s,t) is given by 
and 

" 2{s + t)T 

This average is what comes at finite temperature in place of (E 



(47r)2(s + t)VA 
iN 



The following notations are used [39]: 

H 



(63) 

(64) 
(65) 

(66) 



St 



s + t 



m{s, t) 
A 



1 , u_ 

2 ^+ 

t + sinh(s)e=^'. 



(67) 
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The details on calculations of the PT ( 153|) . (162|) are given in the mentioned 
papers. We note here that the PT is not transversal. Instead, the following 
weaker property holds: pxUxyp'x = 0. The important technical achievement of 
these calculations is the carried out integration by parts, what results in the 
simple explicit expressions for the form factors. 

In the present study, we are interested in two special cases - the Debye, , 
and magnetic, m^agn.y masses for charged gluons. As concerns the Debye mass, 
it coincides with Eqs. fH5l) . fHHj) obtained above for the neutral PT. It obviously 
must hold in the limit of zero background field. However, that it holds for any 
magnetic field strength is to some extend unexpected result. 

Now, we consider the PT in the limit of high temperature, T — oo. It is 
obtained by taking the contribution of the zeroth Matsubara frequency, = 
in the otherwise unchanged expression. This is, of course, the known dimen- 
sional reduction to a theory without temperature in three dimensions. In the 
representation in terms of parametric integrals, the tachyonic projection of the 
polarization tensor takes the form [39] : 



{t\n\t) 




(6^ 



where now f = ll In this expression, we removed the regularization. This is 
possible since the ( 168|) does not contain ultraviolet divergencies. The linear ones 
cancel as in four dimensions and the logarithmic ones do not appear due to the 
dimensionality. 

In the following we focus on the magnetic mass of the charged gluon, 

<agn. = miB,T,p, = 0,p3 -> 0)\t). (69) 

It is given by Eq. f l68|) with = in the tachyonic projection. So it remains 
to treat the tachyonic mode in fl68|) . Because of = this is here simpler than 
in the preceding section. Accounting for the changes dimensionality it can be 
calculated easily and it delivers the imaginary part, 

lm{{t I n I t)T^^) = -1%^ 9'^^^- (70) 



The real part is also much simpler than in the preceding section and it reads 



Aqe-' 2 2 _ 

au — 5 — 2e 

/if_ q smh(g) 



2 TVgB 



-3-26 9' (71) 
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The last line is a result of numerical integration. 

The main result of this calculation is that the ground state projection of the 
PT is proportional to T^/gB allowing for the conclusion that there is a field 
dependent magnetic mass. Further, it is very important that the real part has 
a negative sign. Hence it follows that radiation corrections act to stabilize the 
spectrum of charged gluons at high temperature. The stability parameter is p ~ 1 
as for the neutral gluon case. 

In the electroweak sector of the SM, in the restored phase it is common to 
consider properties of the scalar and gauge fields separately. Here, we follow this 
tradition and say about SU{2) gluons. 

In this and previous section we determined the main screening characteristics 
of gauge fields at high temperature - Debye's and magnetic masses. We observed 
that in the field presence these parameters differ in an essential way from that of 
zero field case. Especially interesting for us here is magnetic mass responsible for 
properties of magnetic fields at high temperature. So, to obtain the results inde- 
pendent of perturbative expansions, in the next two sections we investigate the 
spontaneous magnetization and the neutral gluon magnetic mass by considering 
gauge field theory on the lattice. 

7 Spontaneous vacuum magnetization in lattice 
investigations 

Now, we are going to investigate the spontaneous vacuum magnetization at finite 
temperature in SU{2) gluodynamics on the lattice by using the Monte Carlo 
(MC) simulations. Recent review on the lattice QCD in background fields is 
given in [JO]. In particular, influence of external chromomagnetic fields on the 
deconfinement phase transition was invested in |1T], [12]. It was observed that 
the phase transition temperature, T^, is strongly dependent on the strength of 
the field applied. 

In our problem [20], [21], the main object is a magnetic flux instead of a field 
strength. This is because on the lattice the magnetic field strength is quantized. 
So, it is very difficult to tune the temperature parameter which corresponds to the 
field strength generated inside of the whole lattice. To overcome this difficulty, we 
relate free energy density of the continuous magnetic flux to the effective action 
according to the definition, 

F(^) = S(^,/3)-S(0,/3), (72) 

where S{(p, (3) and ^(O, (3) are the effective lattice actions with and without chro- 
momagnetic field, correspondingly, cp is the field fiux. To detect the spontaneous 
creation of the field it is necessary to show that free energy has a global mini- 
mum at non-zero fiux, (pmin 7^ 0, and its value is negative in the minimum. Let 
us discuss this in more details. 
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As is known, on the lattice free energy is not measured straightforwardly. 
Instead one measures its derivative with respect to the inverse temperature pa- 
rameter P and then free energy is calculated after corresponding integrations. In 
our problem, we, first of all, are searching for the qualitative features of the effect. 
The variable ( !72ll is the corresponding parameter, which has to be measured at 
a given temperature. As we showed above, the field strength generated sponta- 
neously at finite temperature is proportional to coupling constant and therefore 
not strong compared to the squared of the temperature. So, it is expected that 
variable F{ip) takes not large values. Hence, the influence of the flux presence on 
the parameters of the MC lattice calculation in not essential and can be neglected 
in what follows. A detailed investigation of the influence of magnetic fields on 
the MC procedures and parameters is given in [13], [H], where, in particular, it 
is shown that this is really the case. In the course of the standard calculation 
procedures, formula (172]) can be approximated as, 

= c[(5(^, /3) - 5(0, /?)) - (Siif, 0) - 5(0, 0))], 

where c is a positive constant. It is determined by the ratio of the given /3 and 
discrete interval A/3 chosen in the approximation. The terms in the first brackets 
give just cF{ip), and the second brackets is expected to be zero. This is because 
at zero temperature we have a confinement, and the calculated action with and 
without flux is the same. 

In the MC simulations, we use the hypercubic lattice Lt x {Lt < Ls) with 
the hypertorus geometry. and Lg are the temporal and the spatial sizes of the 
lattice, respectively. In the limit of L., — )• oo the temporal size Lt is related to 
physical temperature. 

The standard Wilson action of the SU{2) lattice gauge theory is 



(74) 



= /3E E [l - ^Tr [u^(x)U,(x + afi)\]\{x + au)Vl[ 

where (3 = A/g^ is the lattice coupling constant, (? is a bare gauge coupling, 
U^(x) is the link variable located on the link leaving the lattice site x in the /z-th 
direction. The link variables \J ^{x) are SU{2) matrices decomposed in terms of 
the unity, /, and Pauli Tj, matrices in color space, 

rof., ^ ,^.rpu, - ( + Ul{x) + z^^(x) ^ 



uAx) = im{x) + iTmx) = ( ..a) i ( { . (75) 



Next let us incorporate the external Abelian magnetic field ([2]) into this for- 
malism [20], [l5]. The constant homogeneous external flux (f in the third spatial 
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U„({x.,x„L„x,}+v) = U,({x.,x„0,X3) 

-*>•' — identified — ►<> 



.}+v)e' 



V 



U„(x„x„l 



— identified — ^f- 

,L„x,) = U,(x„x„0,x,)e' 



Figure 6: The plaquette presentation of the twisted boundary conditions. 

direction can be introduced by applying the following twisted boundary condi- 
tions (t.b.c.) |45j: 



The edge links in all directions are identified as usual periodic boundary condi- 
tions except for the links in the second spatial direction for which the additional 
phase ip is added (Fig. [Tj). In the continuum limit, such t.b.c. settle the magnetic 
field with the potential = (0, 0, Bx^, 0) (Ej). The magnetic flux (p is measured 
in angular units and can take continuous values from to 27r. More details on 
the t.b.c. can be found in Ref. [I6j. 

The MC simulations are carried out by means of the heat bath method. The 
lattices 2 x 8^, 2 x 16^ and 4 x 8^ at /3 = 3.0, 5.0 are considered. These values of 
the coupling constant correspond to the deconfinement phase and perturbative 
regime. 

The effective action depends smoothly on the flux ip in the region y9 ~ 0. 
So, the free energy density can be fitted by a quadratic function of (p, 



Uf,{Lt, Xi,X2,X3) 
U^{xo,Ls,X2,X3) 

Ui,{Xo,Xi,X2,Ls) 



Uf,{0,Xi,X2,X3), 
U^{xo, 0,X2,X3), 



e*'^f/^(xo,Xi,0,X3), 

U^{Xo,Xi,X2,0). 



(76) 



These give 




F{ip) = Fmin + Kip - ipminf- 



(77) 
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Figure 7: x^"fit of the free energy density on lattice 2x8^ (left) and 2 x 16'^ (right) 
for ^ = 3.0 (yellow regions describe the ifmin = 0.019l[]:[5}i and ipmin = 0.0069t|j:[]|]P, 
at the 95% confidence level). 



In Eq. f l77|) . there are three unknown parameters, -F^m, b and v'mm- V^mm 
denotes the minimum position of free energy, whereas the -Fmm and h are the free 
energy density at the minimum and the curvature of the free energy function, 
correspondingly. They have been fitted by a standard method. 
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Table 1: The values of the generated fluxes (pmin for different lattices (at the 
95% confidence level). 

The fit results are given in the Table 1. As one can see, (fmin demonstrates 
the 2(T-deviation from zero. The 95% CL. area of the parameters Fmin {b for the 
right figure) and fmin is represented in Fig. [71 The black cross marks the position 
of the maximum-likelihood values of Fmin ip for the right figure) and (^mm- It 
can be seen that the flux is positively determined. The 95% CL. area becomes 
more symmetric with the center at the Fmin-, b and ifmin when the statistics is 
increasing. This also confirms the results of the fitting. 

Thus, we see in this approach that spontaneous vacuum magnetization does 
take place at high temperature. This is in accordance with the results in Sect. 3. 

In Ref. [21] the spontaneous vacuum magnetization with accounting for the 
Ao-condensate at high temperature was investigated. This case has been studied 
in one-loop approximation in [16], [31]. On the lattice, this parameter was in- 
troduced trough the Polyakov loop, as usually. The main result of these studies 
is that the spontaneous vacuum magnetization holds and the vacuum is stable 
for this case. Moreover, with taking into consideration of the condensate, the 
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vacuum stability is increasing. Thus, at high temperature both condensates the 
constant Aq and B(T) are present. 

8 Magnetic mass on the lattice 

In this section we calculate the magnetic mass of the Abelian chromomag- 
netic field by using MC simulations. For that we, following Ref. |16], investigate 
the behavior of the average magnetic fiux penetrating a lattice plaquette oriented 
perpendicular to the magnetic field direction. We introduce the classical mag- 
netic field ([2]) on a lattice as before by applying the twisted boundary conditions. 
Note that in Ref. ^46j the twist of the boundary conditions was applied to in- 
troduce the magnetic fiux of the Dirac monopole. Then, the magnetic mass of 
this non- Abelian magnetic field was measured by investigating the average pla- 
quette values for the twisted and untwisted lattices. The main object of such 
type investigations is the difference (magnetic fiux through a lattice plaquette 
perpendicular to the OZ axis): 

{U untwisted) ~ {Utwisted) = f{f^)Ls), (78) 

which is fitted for each lattice geometry Lt x by different functions /(m, Lg). 

Below we follow this general approach. However, here we measure the mag- 
netic mass of the Abelian field of interest, Eq. (176|l . The temperature is introduced 
as before through a lattice asymmetry in the temporal direction {Lt < Lg). The 
measurements were fulfilled for the value of /3 = 2.6 in the perturbation regime 
for the deconfinement phase. Lattices with = 4 and up to 32 were used. 
To update the lattice, heat-bath algorithm with overrelaxation was used [77] . 
To thermalize the system, up to 6000 MC iterations were used. The plaquette 
average is calculated by averaging up to 10000 working iterations. 

To estimate the behavior of magnetic fields a large amount of simulation data 
must be prepared. Unfortunately, traditional computational resources are lack to 
perform the detail analysis. In our case, we use the General Purpose computa- 
tion on Graphics Processing Units (GPGPU) technology allowing to study large 
lattices on personal computers. GPU programming model implemented here and 
some technical details on MC simulations on AMD/ ATI graphics processing units 
(GPU) are given in Ref. [48j. 

Distinguishing feature of the employed program model is that all necessary 
data for simulations are stored in GPU memory. GPU carries out intermediate 
actions and returns the results to the host program for final data handling and 
output. We avoid any data transfer during the run-time between the host program 
and kernels, to speed-up the execution process. 

To generate the pseudo-random numbers for MC procedure, three different 
pseudo-random number generators are used: RANMAR, RANLUX and XOR128 
|49] . The last one allows to obtain the maximal performance but is not widely 
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used in MC simulations. So, all the results were checked with the slower genera- 
tors RANMAR and RANLUX. 

The GPU-based MC program allows to calculate the difference (|78l) for a wide 
interval of lattice geometries. Also, up to 1000 independent runs for each lattice 
size were performed in order to decrease the dispersion of the obtained values 
f{m,Ls). The whole set of simulation data for different lattice geometries were 
fitted with the several functions which correspond to the different behavior of 
magnetic flux. 

The results of fitting (fitting function, the values of the x^-function corre- 
sponding to the 95% confidence level and the obtained magnetic mass m) are 
shown in Table 2. The function ^ corresponds to the magnetic flux tube for- 
mation (r is the lattice size in the X and Y directions). The total magnetic 
flux through the lattice is conserved in this case. The function ^ describes the 
Coulomb-like behavior and the function ^exp(— m^r^) is signaling the genera- 
tion of the magnetic mass m [IB]. The functions ^exp(— mr), ^ can be related 
to the increase of the field strength with a temperature increase. This is because 
the total magnetic flux through the lattice is growing faster than in the case of 
the magnetic flux tube formation. 

As it follows from Table 2, the best fit function is — exp(— mr) with a small 
value of the magnetic mass m = 1.25 ■ 10~^. The value of function in this 
case is very close to the m = situation and statistically these cases are indis- 
tinguishable. Really, the statistical errors are larger than the fitted value of m. 
Thus, from the carried out analysis we can conclude that the neutral component 
of the gluon field is not screened at high temperature like usual magnetic field. 
This result is in agreement with that of section 4 obtained in perturbation theory. 
Note that due to a large amount of data we have guarantied that the absolute 
value of errors is of 10^ — 10^ times smaller than the value of the corresponding 
quantity. 
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Abelian field 


Fit function 




C 


■parameter 


C exp(— mr) 


901.8 


0.063 


m = (2.44^0:06) . 10-2 


C exp(— m^r^) 
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0.035 


m = {i.hitri) ■ 10-2 


C/r 
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0.911 




C/r exp{—mr) 
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m = (1.25t^4) • 10"^ 


C /r exp{—m?r'^) 
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^2 = (2.4l|?51.2) . 10-10 


cy 


31400 


28.13 




Cjr^ exp(-m2r2) 


7550 


18.26 


m? = -3.3 • 10-5 




159500 


248.9 




C/r^ exp(— m^r^) 


161000 


10.0 


m = 0.0 



Table 2: Fit results for magnetic mass of Abelian magnetic field. 

Interesting additional arguments in favor of spontaneous vacuum magnetiza- 
tion follow from the above measurements fulfilled. We observed that for the fitting 
function /(m, L^) = ^ corresponding to the magnetic flux tube formation the 
value is very large and entirely inconsistent with the data. But in the geometry 
of mcasm^cmcnts it describes the conservation of the magnetic flux introduced 
by the twist of the boundary conditions. The best flt functions ^,^exp(— mr) 
with very small (actually, zero) m arc signaling an increase of the mean magnetic 
fleld strength penetrating the plaquettc perpendicular to the fleld direction. As 
a result, the flux through the whole {X — Y) plain should increase. The only 
natural explanation is the spontaneous generation of the fleld inside the volume 
of the lattice. 

9 Magnetic field strength at Tg^^ 

In previous sections, we have shown that the spontaneous vacuum magnetiza- 
tion is realized at high temperature and determined the screening parameters 
mD{B,T) and rrimagn.- Now, we are going to estimate in the SM the magnetic fleld 
strength at the temperatures close and higher the EWPT temperature, T > T^w 
As it was shown in Sect. 4, in the restored phase the hypermagnetic fleld 
h^j, — {) and the complete weak-isospin chromomagnetic fleld A^^^ is generated. 
It is unscreened because its magnetic mass is zero. Thus, the fleld is a long- 
range one. It provides the coherence length \b(T) to be sufficiently large. In 
fact, the fleld is a constant and occupies all a horizon scale at given temperature, 
Ab(T) ~ Rh(t)- Chromo (magnetic) flelds of different types (color 5^/(3), and 
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others) can be spontaneously generated at high temperatures. This property 
could be of great importance for cosmology. 

Now, we continue with describing a general field behavior related with the 
EWPT. In the restored phase, a scalar field condensate = and the con- 
stituent of the weak isospin field corresponding to the magnetic one is given by 
the expression 

B{T) = sin 9UT)B^^\T), (79) 

where B^^\T) = B is the strength of the field generated spontaneously. After the 
phase transition, the scalar condensate (p ^ and the field is partially screened. 

To estimate the field strength B{T) in the restored phase at the EWPT tem- 
perature the total EP must be used. This can be best done numerically. To 
explain the procedure, we consider here the part of this potential accounting for 
the one-loop ly-boson contributions. The high temperature expansion for the 
EP coming from charged vector fields is given in Eq. ([71). Assuming stability of 
the vacuum state, we calculate the value of chromomagnetic weak isospin field 
spontaneously generated at high temperature from Eqs. ([7]) and ([8]): 



1^ 

16 7r2(i + A^logI; 



gB{T) = — ^ — —. (80) 



Here, r is a reference temperature parameter. This expression (and the com- 
plete one accounting for all the contributions) gives the field strength at any 
temperature T > Tew Such type formulas can be obtained for different models 
of particles. 

Now, before educing a specific value of it, we describe how to relate this 
expression with the present day intergalactic magnetic field Bq. We first relate 
the expression flHOj) with an magnetic field after symmetry breaking, and then 
take into account the scales of fields. 

2 

Let us introduce the standard parameters and definitions, |^ = a^, a = 

as sin 6*^, = ay and tan^ Oy^iT) = , where a is a fine structure constant. 
To find the temperature dependence of the Weinberg angle, the behavior of the 
hypercharge coupling g' on the temperature has to be computed. From Eq. it 
follows that this behavior is not trivial. The logarithmic temperature-dependent 
term is negative. But, as it is well known, in the asymptotically free models this 
sign must be changed to a positive value due to contributions of other fields. 
This particular value is model dependent. We will not calculate it in the present 
paper. Instead, for a rough estimate, we substitute the zero temperature number: 
sin^e^iT) = sin^e^iO) = 0.23. 

For the given temperature of the EWPT, Tew, the magnetic field is 

B{Tew) = Bo^ = sinew{Tew)B{Tew). (81) 
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This relation is the consequence of the assumption that for the field sponta- 
neously generated at high temperature in the early Universe the magnetic flux 
conservation holds after the EWPT. It means that the field is "frozen" in plasma 
at large scales and the magnetic turbulence processes do not affect this behav- 
ior. Although this is most simple assumption, it requires a detailed discussion 
within the results obtained recently in magnetic hydrodynamics (see [50], [51] 
and references therein). 

Assuming T^^ = lOOGeV = lO^^eV and Tq = 2.7K = 2.3267 • 10~^eV, we 
obtain 

S(Te^) ~ 1.85 ■ lO^^G. (82) 

This value can be considered as an estimate of the magnetic field strength at 
the EWPT, if the standard model serves as a basic theory describing this phase 
transition. Hence, for the value oi X = log ^ we have the equation 

° " 2 7rV2sin2^^^(l + ^-5%-X)2- ^^^^ 

Since all the values are known, logr can be computed. After that the field 
strengths at different higher temperatures can be found. In fact, the main point in 
obtaining of these results is the assumption on the magnetic flux conservation and 
frozen in it in the plasma. An information on a particular model is implemented 
in the factor smOwiTew) in Eq. (IHT]) . 

Of course, our estimate is a rough one because of having ignored the temper- 
ature dependence of the Weinberg angle. To guess the value of the parameter 
r we take the field strength Bq ~ lO^^G, often used in the literature (see, for 
example, [S2])- In this case, from Eq. ( 15^ we obtain r ~ 300eV^. For the present 
day value Bq ~ 10~^^G this parameter is much smaller. 

Let us compare now the value of the field strength fl82|) with the one calculated 
directly from the EP for the SM in Ref. [53]. From Fig. 1 and Tab. 1 of that 
paper, we find 

B'^{T,^) ~ 102°G, (84) 

what is much larger than the value ( 182|) and just corresponds to the value of 
the present comoving field strength Bq ~ 10~^G. This value was considered in 
numerous investigations as an upper bound on the intergalactic magnetic field 
strength in the Universe, before the recent discoveries [in]-[I2]- Note that the 
value of B(Tew) calculated with expression (150]) is close to the estimate 
Let us stress again that the field strength at higher temperatures will depend on 
the particular model extending the standard one. Spontaneous vacuum magne- 
tization in the minimal supersymmetric standard model has been investigated in 
Ref. [30], and the field strength generated in this model is smaller as compared 
to the situation here considered. Also, Pollock [52] has investigated this problem 
for the case of the Planck era, where magnetic fields of the order B ~ lO^^G have 
been estimated. 
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An important consequence of Eq. (IHO!) is that at high temperatures T > Tg^ 
the magnetic flux conservation does not hold. This follows from the logarithm 
term standing in the denominator. Therefore, at a given temperature T a specific 
magnetic field strength is generated due to vacuum polarization. This property 
is very important for what follows. It is a key point for the scenario which we 
are going to investigate. We assume for the following that intergalactic magnetic 
fields have been spontaneously generated in the hot Universe. 

Note also that magnetohydrodynamical processes in the early Universe were 
investigated in numerous publications (see for references the review papers [5^ . 
[T5]). Here, we mention only the important points for our consideration. The 
"frozen in" conditions can be realized for magnetic fields having the scales of 
largest turbulence eddies. After a free decay stage of the magnetized plasma 
evolution, the field can be considered as non-affected by turbulence [51]. In 
connection with these results, it is clear that the magnetic fields generated at the 
EWPT are not sufficient as such to produce a long-range correlated fields and 
some additional processes must be included. This is because even the field having 
the scale of the Hubble radius at temperature Tg^ ~ 100 GeV is correlated at the 
comoving scale /q ~ 10~^ Mpc [51]. However, the scale of intergalactic magnetic 
fields was determined to be of the order ~ 1 Mpc [10] - [12]. The fields 
generated at the EWPT can serve as seeds for long-range magnetic fields inside 
galaxies. 

10 Magnetic field scale 

Now we discuss the scale of the field generated in the restored phase [55]. This 
is a key point in relating expressions like Eqs. (180|1 or ( 18^ with Bq. In our con- 
sideration, the "frozen in" condition was used. So, let us discuss its apphcability 
in more detail. Remind, if one assumes that after the EWPT the constant field 
B{Tew) was frozen in the plasma at the Hubble scale, Rh{T(.w)i then its comov- 
ing coherence scale at present has to be \b{Tq) = 6 • 10""^ pc pT]. This is much 
smaller than necessary. 

We consider two, in fact related, ways to overcome such a difficulty [55]. The 
first is to take into consideration the reheating stage of the universe evolution. 
According to the concepts of modern cosmology [56], this stage has existed just 
after infiation and is related with the latter causal stage. Just due to causality, 
the temperature in the Universe after this stage is the same, in all domains of 
space, which could even be uncorrelated in later moments of time. Hence, at 
a given high temperature, T, a magnetic field generated due to vacuum polar- 
ization has the same strength B{T) everywhere in the Universe. Formally, the 
field strength could have different directions, in either external or internal spaces. 
Different kind of (chromo)magnetic fields, of the type as in Eq. ([2]), can be spon- 
taneously generated. Their nature depends on a particular model considered and 
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is therefore unknown, so far. But this is not essential for our consideration, here. 
The magnetic fields coherent on huge scales have been present in the early Uni- 
verse. The origin of this coherence is ensured by the properties of the solution 
to the field equations (Eq. ([2])) and by causality at the inflationary epoch. The 
scales of the coherent field domains could be estimated on grounds of the gauge 
invariance. This idea, due to Feynman, was put in force in gluodynamics with 
the goal to determine possible magnetic vacuum structures [57]. Namely, to find 
a gauge invariant vacuum, on the basis of gauge non-invariant solutions (such as 
Eq- (E])), one can consider a domain structure ensuring gauge invariance when a 
corresponding boundary is going around. At T 7^ 0, this point requires further 
investigation. 

Most of fields generated in the early Universe decouple and are screened at 
some energy (temperature) scales, when the corresponding scalar condensates 
have broken the background symmetries. So, the only unbroken symmetry at the 
EWPT remains the SU{2)e^ x f/(l) one. After the EWPT, when spontaneous 
magnetization stops, this field cannot be included in turbulent processes gener- 
ated by the transition. This is because the scale of the field, Eq. ([2]), is already 
much larger than any largest eddy of turbulence. As it is usually believed, the 
size of a typical eddy is estimated as the inverse mass of the particles appearing 
after the transition [58]. Thus, the field evolves in accordance with the metric 
expansion and is implemented in a hot plasma, thus fulfilling the magnetic flux 
conservation law. And it finally results in the present day intergalactic mag- 
netic fields which could be correlated on ~ 1 Mpc scales. Note that an essential 
information on the processes that take place after the EWPT, obtained in the 
framework of magnetohydrodynamics, is given in Refs. [51] and [5U] . 

Another possible scenario is based on the stochastic processes considered al- 
ready by Hogan [60] in connection with the magnetic fields generated at first-order 
EWPT. A possible mechanism of field generation in that case was proposed by 
Vachaspati [58]. In the former paper, it was pointed out that magnetic fields 
correlated on large scales can be produced not only through causal processes 
but also by a stochastic random walk mechanism, if the magnetic lines gener- 
ated in some domain of space "forget" about their origin. The field strength 
developed on large scales by this process (due to "straightening" of entangled 
magnetic fields) can be estimated as ~ B/\fN, where counts the number 
of domains, with the field 5 of a given size, crossed by a magnetic line. The 
correlation length in this case can be much larger than the Rh{T). It can 
be estimated as Xb{T) ~ NRh{T). In Ref. [60], it has been also noticed that 
this mechanism is not applicable to the early Universe, the reason being because 
magnetic lines do not penetrate freely though the plasma. This is really the case, 
if the properties of the plasma as taken into account. However, this is not the case 
if spontaneous vacuum magnetization occurs. Actually, at a given temperature, 
each uncorrelated domain of space having a Hubble radius Rh (T) is filled up with 
a constant magnetic field B(T), described by the potential ([2]). Its orientation 
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in both external and internal spaces is arbitrary. Hence, a stochastic behavior 
of the field lines and the appearance of magnetic fields having large correlation 
lengths Xb{T) > Rh(T) are expected. After the EWPT, these fields evolve as in 
the previous case. 

Note that, in both scenarios, all the fields generated at the infiation epoch [6] 
are washed out by the vacuum polarization and leave no remnants at present. The 
reheating stage becomes more important. In closing this section, we also remind 
that the long-range nature of the Abelian spontaneously generated magnetic fields 
is ensured by their zero magnetic mass, what renders these fields unscreened, as 
is the case for usual U{1) magnetic fields. The difference, however, is essential 
because the former fields appear due to vacuum polarization and the latter ones 
need currents to be produced. 

11 Discussion and conclusions 

We here summarize our main results. The spontaneous vacuum magnetization 
process eliminates the magnetic fiux conservation principle at high temperature. 
Vacuum polarization is responsible for the value of the field strength B{T) at 
each temperature and serves as its source. Physically, the magnetization is a 
consequence of the large magnetic moment of charged gauge bosons (gyromag- 
netic ratio 7 = 2). This is one of the basic properties of matter. In fact, at high 
temperature, due to spin interaction of the virtual pairs, heat is converted into 
an ordered coherent vacuum state. This effect was not taken into consideration 
in previous studies of the early Universe. 

We have shown here that, at finite temperature and after symmetry breaking, 
a scalar field condensate suppresses the magnetization. Hence, it follows that the 
actual nature of the model extending the standard one is not that essential at 
sufficiently low temperatures, when the decoupling of the heavy gauge fields has 
occurred. From this one can conclude, in particular, that the vacuum polarization 
"washes out " the relics of the magnetic fields generated at very high temperature 
or at the infiation stage. This is because different kinds of magnetic fields existed 
at high temperatures were screened with temperature lowering, and the magnetic 
flux conservation does not hold at T > T^w After the EWPT, when the solution 
7^ is realized, the spontaneous vacuum magnetization stops and the field 
evolves further according to the magnetic fiux conservation law. The correlation 
length of the field is very large. At the EWPT, the field can be considered as 
constant magnetic field having the strength B{Tew) ~ 10^^ G. This value is not 
large and the field not essentially infiuences various processes happening after the 
transition. This is because, according to the estimates in Ref. [58], the stochastic 
magnetic fields generated in a first order phase transition are much stronger, 
B ~ 10^° — 10^'^ G. The scale of these fields is of the order of inverse particle 
mass, A ~ l/fi'^- 
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Recall that after reheating, the temperature in the hot Universe is the same 
everywhere (even in the domains of space which are uncorrelated in later periods 
of time). Analogous statement concerns the values of B{T) which are strictly 
related with the temperature due to the vacuum polarization. Most magnetic 
fields existed at high temperatures were screened with the temperature lowering 
when the corresponding symmetries have been spontaneously broken. The only 
magnetic field at the temperatures close Ten; is the SU{2)ew component described 
by the potential ([2]). 

The present value of the intergalactic magnetic field strength, Bo, we related 
with the field in the restored phase. Because of the zero magnetic mass for 
Abelian magnetic fields, there is no problem for the generation of fields having 
large coherence scales. In our estimates, we have assumed that, basically, the 
field is of the order of the horizon scale, Ab(T) ~ Rh{t)- This seems reasonable 
because, at a given temperature, the field B(T) = const, generated due to vacuum 
polarization, occupies all space. Then coherence on scales exceeding that of the 
horizon can be produced. Such large scale fields are not influenced by turbulent 
processes happening after the EWPT. They are frozen in the plasma and evolve 
according to the magnetic flux conservation law. In this scenario, a large scale 
domain structure is also permissible, what requires additional consideration. 

As we have found above, the field strengths at the EWPT temperature, esti- 
mated with account to the present-day value of the intergalactic magnetic field 
strength, Bq ~ 10^^^ G, Eq. (!82|) . or either directly from the vacuum magneti- 
zation in the standard model, Eq. (18^ . differ in six orders of magnitude. This 
huge deviation can be explained by the different scales of the fields considered. 
Let us check this possibility by using the second of the scenarios proposed in 
the previous section, for large scale field generation. We used the usual relation 
between the scale factor and the temperature, 

(85) 

taken at the EWPT epoch, and the present-day parameters. Tew = 100 GeV = 
10^^ eV, To = 2.3267 ■ 10"^ eV. If one now assumes that Ab(T) ~ a(T), then 
from (!85l) it follows that Xb(To) = 6 ■ 10~^ ps. On the other hand, if one takes 
A_b(To) = 1 Mpc, the value XB(Tew) = 2.33 ■ 10~^^ Mpc is obtained. At the same 
time, the horizon size is a(Teu,) = 1.27 ■ 10"^'^ Mpc, thus, XB(Tew) » aiTew)- 
Now, following an idea of Hogan [60], we relate the size of the correlated field 
with the random walk process. At Tew, we have XsiTew) = Na{Tew), hence, we 
get roughly -y/iV = 3 ■ 10^, and for the field strength "straightened" on the A^- 
domain scale, -Bjv ~ ^^^^ |60] . Therefore, accounting for the field strength value 
calculated for the standard model, Eq. fl8^ . we obtain BisiTew) ~ 3 ■ 10^^ G (the 
subscript in Bi^ means "large scale"). This value is close the value BigiTew) ~ 
2 ■ 10^^ G estimated in Eq. ( I82p . The remaining discrepancy can be explained in 
two ways. First, and obviously, as due to the roughness of our estimate. Second, 
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and more radically, by the necessity of substituting the standard model with 
another one. The latter point will be discussed in more detail below. 

Let us note a number of properties of the field under consideration, and com- 
pare them with the ones usually discussed in applications of magnetohydrody- 
namics to the early Universe. Here, we follow Refs. [51j and [59] which are close 
to our analysis. First, let us mention that the field energy density, pb = is 
proportional to g^T'^, what is much smaller than the radiation energy density, 
~ T'^. Thus, the BBN condition (see [51]), pb/ Prad. « 1, is fulfilled. Second, as 
numerical simulations show [Slj, the turbulent process in the early Universe with 
magnetic field included is slower than in the laboratory. Turbulence can include 
a large scale field at the level of the largest eddies. For large scale fields, the 
free decay stage is important. At this epoch, which is strongly dependent on the 
initial conditions [59], turbulence is significantly decreasing, and after this very 
brief stage the field is not affected by turbulence any more. It is just frozen in the 
plasma. As we have shown, the spontaneous vacuum magnetization is stopped 
when the first-order phase transition ends. That is, just after the free decay stage. 
Fields of this type cannot be influenced by turbulence and thus it is reasonable 
to believe that after the EWPT the field evolves according to Eq. flHTl) . Note also 
that these fields are non-helical ones. 

Now, we compare our results with those of Ref . , where spontaneous vac- 
uum magnetization at high temperature was applied to estimate the field strength 
at the Planck era. In that paper, to estimate the value B{T) the heterotic super- 
string theory x was considered as the basic model. At the Planck era, the 
magnetic field strength has been estimated to be of order ~ lO^^G. In contrast to 
our considerations, it was assumed there that the magnetic field approximately 
scales as -B ~ T^. That is, vacuum magnetization was taken into account only at 
the very first moments of the evolution. Further, recent results deriving a zero 
magnetic mass for the Abelian chromomagnetic fields also change the picture 
of the magnetized early Universe substantially. According to this, the created 
magnetic fields existed on the horizon scales. They were switched of at some 
mass scales. The processes of decoupling were also not taken into consideration 
in Ref. [52]. Thus, it is impossible there to relate the magnetic field Bq with the 
magnetic fields generated at high temperatures. 

In the used approach, knowing the particular properties of the extended 
model, it is possible to estimate the field strengths at any temperature. This 
can be done for different schemes of spontaneous symmetry breaking (restora- 
tion) by taking into account the fact that, after the decoupling of some massive 
gauge fields, the corresponding magnetic fields are screened. Thus, the higher the 
temperature, the larger the number of strong long-range magnetic fields of dif- 
ferent types that will be generated in the early Universe. To estimate their field 
strengths at higher temperatures, one has to take into consideration a number of 
features proper to the standard model and its particular extension at play. First, 
we note that quarks possess both electric and color charges. Therefore, there is 
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a mixing between the color and usual magnetic fields owing to the quark loops. 
Second, there are peculiarities related with the particular content of the extended 
models. For example, in the Two-Higgs-Doublet standard model the contribution 
~ {gBY^'^T in Eq. (JTj) is exactly canceled by the corresponding term in Eq. (jH]), 
because of the four charged scalar fields entering the model. They interact with 
gauge fields with the same coupling constant. However, in this model the dou- 
blets interact differently with fermions. This changes the effective couplings of 
the doublets with the gauge fields and results in non-complete cancelations. As 
a result, a suppression of the spontaneously created magnetic field is expected 
in this model. In principle, one should be able to explain, in this way, the dis- 
crepancy in the field strengths as discussed above. Note also that in the minimal 
supersymmetric standard model the field B{T) generated at high temperature is 
smaller then in the standard model [30]. So, formally, it can also be considered 
as a candidate for explanation of intergalactic magnetic fields. However, this case 
is not well studied, now. There can be other peculiarities which may influence 
the high temperature phase of the Universe. 

As conclusion we note that the spontaneous vacuum magnetization at high 
temperature could find other interesting applications in cosmology. Another phe- 
nomenons of interest may be a deconfinement phase transition with accounting 
for both magnetic B{T) and color chromomagnetic fields Bc{T) and a chiral 
magnetic effect. These require further investigations. 
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